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Abstract 



We investigate whether the (planar, two complex scalar) dilatation operator of M = 4 
gauge theory can be, perturbatively and, perhaps, non-perturbatively, described by an 
integrable long range spin chain with elliptic exchange interaction. Such a chain was 
introduced some time ago by Inozemtsev. In the limit of sufficiently "long" operators a 
Bethe ansatz exists, which we apply at the perturbative two- and three-loop level. Spec- 
tacular agreement is found with spinning string predictions of Frolov and Tseytlin for 
the two-loop energies of certain large charge operators. However, we then go on to show 
that the agreement between perturbative gauge theory and semi-classical string theory 
begins to break down, in a subtle fashion, at the three-loop level. This corroborates 
a recently found disagreement between three-loop gauge theory and near plane-wave 
string theory results, and quantitatively explains a previously obtained puzzling devia- 
tion between the string proposal and a numerical extrapolation of finite size three-loop 
anomalous dimensions. At four loops and beyond, we find that the Inozemtsev chain 
exhibits a generic breakdown of perturbative BMN scaling. However, our proposal is not 
necessarily limited to perturbation theory, and one would hope that the string theory 
results can be recovered from the Inozemtsev chain at strong 't Hooft coupling. 



1 Introduction and conclusions 



Nearly two years ago a fresh approach pP to uncover the dynamics of the AdS/CFT 
correspondence [2j was proposed. Progress in this direction is clearly required if we 
are to understand the gauge theory implications of the correspondence for quantum 
gravity and string theory, as well as the string theory implications of AdS/CFT for 
gauge theory. Berenstein, Maldacena and Nastase (BMN) [T] argued that the string 
spectrum on AdS^ x S 5 in a linearized (plane wave) limit, which is known when 
expressed in planar Af = 4 gauge theory language, leads to a prediction for the all- 
loop exact anomalous dimensions of certain operators containing one large R-symmetry 
charge J. These predictions take the form of explicit functions in the BMN coupling 
constant 

A' = ± (1) 

which are analytic around A' = 0, with a finite radius of convergence. One can therefore 
expand them in integer power series in A', and thus, so it seems, A. This led BMN 
to suggest that the predictions could be checked, order by order, in perturbative planar 
M = 4 gauge theory. The claim is that thereby one avoids a notorious AdS/CFT 
difficulty: the regime where AdS string theory is calculable corresponds to A large, while 
gauge theoretic perturbation theory obviously assumes A to be small. We will have more 
to say about this below. 

The BMN proposal triggered a large number of interesting research papers, most of 
which cannot be discussed here (see the thorough reviews [1]). In particular, powerful 
techniques for the efficient computation of one-loop anomalous dimensions of Af = 4 
composite operators containing a large number of scalar fields were developed jS]. These 
activities led to the highly important insight of Minahan and Zarembo [6J that the 
planar one-loop anomalous dimensions of these operators can be found by diagonalizing 
the Hamiltonian of an equivalent integrable spin chain by the Bethe ansatz method. In 
the case of just two complex scalars the spin chain is extremely simple, the Heisenberg 
XXX model. This is the "harmonic oscillator" of condensed matter theory. Similar spin 
chain techniques had been previously developed in the QCD context for a different class 
of operators, rather related to the space-time symmetries instead of the R-symmetry of 
the scalar operators [2] . A unified one-loop treatment of all conformal operators in the 
Af = 4 theory was developed in [H], and shown to yield an integrable PSU(2,2|4) super 
spin chain as well as its corresponding Bethe ansatz. 

The discovery of integrability in perturbative Af = 4 theory is, potentially, of great 
importance. In a parallel development integrable structures were also pointed out on 
the string side of the correspondence [TO]. Since the structures in gauge and string 
theory appear in very different regimes, the exciting prospect arises that the planar 
AdS/CFT system might be completely and non-perturbatively integrable. If true, this 
might allow to "exactly solve" planar Af = 4 gauge theory, and thus free IIB string 
theory on AdS$ x S 5 . 

One way of looking at the BMN approach is to consider it a semi-classical limit jTT] 
with an "artificial" large parameter, the charge J. In fact, a beautiful generalization of 
the BMN limit was proposed by Frolov and Tseytlin ^2] (see also and the closely 
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related earlier work |14j). They consider rotating string solutions of the classical string 
equations of motion with several large charges Ji, J2 ■ ■ ■■ They then show that, firstly, 
quantum sigma-model corrections are suppressed by powers of the inverse total charge 
J = J\ + J 2 + . . ., and that, secondly, the obtained expressions for the string energies 
again expand in integer powers of the BMN coupling constant eq.Q. The latter fact led 
Frolov and Tseytlin to propose that, just as in the BMN case, the conjectured dual gauge 
theory operators should possess anomalous dimensions whose perturbative expansions 
match the just mentioned string energy expansions. An added benefit of this suggestion 
is that the symmetry charges Ji allow for a simple identification of those conformal 
M = 4 scaling operators which are natural candidates for the gauge theory analogs of 
the string theory states. The proposal was further elaborated on a large number of 
explicit examples [Eh [EI- in all cases the classical problem can be exactly solved, the 
underlying reason being, once again, integrability ^7] (classical, however, in this case). 

The final (all-orders in A') expressions for the classical string energies are rather 
non-trivial functions of the charge ratios Ji/J, as they are obtained by solving a non- 
linear system of equations. This is even true for the leading O(X') term, which should 
correspond to a one-loop anomalous dimension in gauge theory. By applying the above 
mentioned Bethe ansatz to the relevant operators, it was indeed shown in a number of 
non-trivial examples that the leading string predictions can be reproduced by one-loop 
J\f = 4 gauge theory [E] , [E] , [2D] • The relevant quantum spin chain is integrable since 
it possesses an infinite number of mutually commuting charges. As was recently shown 
in [21], these may be reproduced from the Bethe ansatz solutions as well, and precisely 
match the corresponding infinite tower of string sigma-model charges to leading order in 
A'. For a recent, up-to-date review of "spinning strings", see [22]- 

The last result [21] would seem to be a near-proof, assuming the correctness of 
AdS/CFT, that (1) J\f = 4 is integrable to all orders in A', and that (2) the BMN 
and Frolov-Tseytlin (FT) proposals should indeed be valid to all orders in perturbation 
theory. The charges certainly are commuting for all values of A' on the string side, 
and it would appear to be very difficult to imagine that there exist two inequivalent 
deformations of the matching leading order O(X') commuting charges. Some caution is 
nevertheless warranted. Curious, unexpected and unaccounted structural agreements be- 
tween small A gauge theory and large A string theory results have previously appeared in 
the context of AdS/CFT, see in particular [23] • One should definitely push perturbative 
M = 4 to higher orders and see whether the agreement with BMN/FT persists. 

The BMN prediction, originally checked in [T| at the one-loop level, was successfully 
tested in [23], [22] at two loops. It is usually stated in the recent literature that in [26 j 
an all-orders proof was given. It appears to be difficult to rigorously justify some of the 
details of the proof, and we remain agnostic with regards to its validity. At any rate, 
while inspiring, the methodology of 25] is not of a constructive nature and we are unsure 
whether it does not, just like string theory, implicitly assume large A. 

In [2Z] a program was begun to derive the M = 4 dilatation operator beyond the one- 
loop level. The two-loop correction for two complex scalars was found, and interpreted 
as a next-nearest neighbor deformation of the nearest neighbor Heisenberg spin chain. 
Excitingly, the result showed that the two-loop terms do not lift certain spectral planar 
one-loop degeneracies, which could be traced to the integrability property of [§]. This 
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led to the conjecture j2Zj that the exact dilatation operator might correspond to an 
integrable long-range spin chain. Assuming integrability, and BMN scaling, the three- 
loop correction was derived. Using its explicit form, a numerical two- and three-loop 
gauge theory estimate for the FT prediction for the "folded string" jTHj (see chapter 3 
below) was obtained in ^Hj- It showed excellent agreement at two loops, but a strange 
17% deviation at three loops. As a second sign of potential three-loop trouble, a string 
calculation of the leading 1/ J correction to the BMN anomalous dimension formula was 
performed in [2E] and, while beautifully agreeing at two loops, it cfoagreed at three loops 
with [27]. 

Important further information comes from recent work by Beisert j2H]- The method 
goes beyond [2TJ and uses, apart from field theoretic structural properties, symmetry 
arguments in order to constrain the three-loop dilatation operator (actually, for a larger 
class of fields than just two complex scalar s). For our present purposes (we stick to the 
latter case in the present paper), the upshot of his study is as follows: (1) Perturbative 
integrability (in the sense of |2?| ) extends to at least three loops. (2) The two-loop 
dilatation operator is fixed up to one constant. This constant may be fitted to the 
known two-loop anomalous dimension of the Konishi field. This way one independently 
verifies the two-loop BMN prediction, in agreement with [21], [2H]- (3) The three- loop 
dilatation operator is determined up to two unknown constants. These could be found 
if we knew two (independent) M = 4 anomalous dimensions (such as Konishi, plus 
one further field). However, no three-loop Af = 4 anomalous dimensions are rigorously 
known to date 1 . Therefore, the three-loop dilatation operator is not yet rigorously known. 
However, if we impose qualitative BMN scaling (i.e. if we assume in accordance with 
eq.Q that the 0(X 3 ) anomalous dimensions of finite J BMN operators [3U] scale like 
~ J~ 6 ), both constants are fixed and the infinite J BMN result is also quantitatively 
reproduced 2 . 

In an interesting parallel development, it was established that a system closely related 
to full-fledged Af = 4 gauge theory, plane-wave matrix theory, also exhibits integrability 
to at least three loops j^I]. This is a quantum mechanical system, and this result was 
rigorously derived without further assumptions. 

Summarizing the last paragraphs, integrability seems to be a very stable concept in 
Af = 4 gauge theory even beyond one-loop; with further, indirect evidence coming from 
the string side of the AdS/CFT correspondence [TQJ , [TJJ , [21] • It is therefore important to 
gain deeper understanding of the integrable long range spin chain potentially describing 
planar Af = 4 gauge theory non-perturbatively. Such a spin chain, in contrast to the 
Heisenberg model, should contain an extra parameter, related to the Yang-Mills coupling 
constant. There exists one well-known integrable long-range spin chain, namely the 
Haldane-Shastry chain [22] • However, there the parameter corresponds to the length of 
the chain 3 . In particular, it is not possible to recover the one-loop Heisenberg model from 
Haldane-Shastry. There also exists a much less known integrable spin chain that precisely 

1 It would clearly be of great importance to perform a reliable field theoretic computation of the 
planar three-loop anomalous dimensions of two independent fields in J\f = 4 gauge theory. 

2 However, as we just mentioned, the 1/J corrections are not properly reproduced 28 . 

3 This is reminiscent of the BMN coupling constant eq.QJ, which also depends on the length of the 
spin chain. From our point of view, this is not an accident. 
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furnishes such an extra deformation parameter, the Inozemtsev long range spin chain [33J. 
Interestingly, it encompasses both the Haldane-Shastry chain and the Heisenberg chain. 
It is the purpose of this paper to investigate whether it might serve as a candidate for 
describing the sought all-loop dilatation operator (for two complex scalars). We shall 
find that this is indeed consistent with all currently known solid information. However, 
proving or disproving it will require more work. 

The paper is organized as follows. In chapter two we describe the Inozemtsev spin 
chain, and show that it is indeed capable of emulating the M = 4 dilatation operator up 
to at least three loops. It does not fix the yet to be determined two constants of the three- 
loop terms, as discussed above. This is not surprising, as they do not affect integrability 
[2*H] . We also describe a Bethe ansatz for this quantum chain, which is due to Inozemtsev 
as well (see [HI], and references therein). This ansatz (unlike the full Inozemtsev spin 
chain) is limited to operators which are "long" , which means that the interaction range 
is kept smaller than the length of the chain. We employ it to re-derive all known results 
|2f] for finite J BMN operators. The Inozemtsev-Bethe ansatz is particularly interesting 
since it may, in principle, be used to derive non-perturbative results. In chapter three we 
apply the ansatz to the case of spinning strings. We are able to derive some of the highly 
involved two-loop expressions of [T2"l ITol IT7] for folded and circular strings, in agreement 
with our previous numerical study for a particular case [T^]. However, at three loops 
we find that gauge theory yields a similar, but different result as compared to the FT 
proposal. On the positive side, we very accurately reproduce the previous numerical 
estimate for a special case [THj . 

Is this really bad news? Not necessarily. Of course, the perspicacious reader might 
suspect that the three-loop dilatation operator as conjectured in [27] and derived in j2H] 
is simply wrong. In the absence of an independent, rigorous field theory computation, 
this remains a theoretical possibility; however, we currently do not see which mistake 
could have been made. On the other hand, it might be the case that the BMN and FT 
results are only valid at large A (the one- and two-loop agreement might be due to some 
yet-to-be understood planar weak-strong coupling duality symmetry); the expansion is, 
after all, in A', not A. This issue came up already in certain puzzles related to the 
BMN limit [35 3 . In this context it is interesting that the Inozemtsev chain exhibits a 
generic breakdown of perturbative BMN scaling at four loops 4 . However, it is possible 
to recover BMN scaling non-perturbatively, i.e. at strong coupling (see chapter 2). A 
further, more troubling explanation for the deviation we are finding would be a large iV 
phase transition when we go from weak to strong coupling [36]. 

In this paper we have shown how to extend the Bethe ansatz for planar one-loop 

4 We note in passing that we suspect the Inozemtsev model to be equally capable of describing the 
plane- wave matrix model mentioned above, which has also been shown to be "perturbatively" integrable 
to at least three loops There BMN scaling may be recovered, up to three loops, by a redefinition of 
the mass parameter of the model. If our Inozemtsev ansatz is correct we suspect that such a redefinition 
is no longer possible at the four-loop level (or else, integrability of the model should break down). This 
might be a more feasible check of our Inozemtsev ansatz as opposed to considering four- loop perturbation 
theory in full-fledged J\f = 4 field theory. In turn, if the plane wave mechanics stays integrable at four 
loops while still being consistent, after redefinition, with BMN scaling, it would strongly suggest that 
there exist further integrable long range spin chains differing from Inozemtsev. 
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J\f = 4 gauge theory to higher loops. Our procedure is completely sound up to three 
loops. Of course we are far from proving that our proposal is valid beyond three loops, 
let alone non-perturbatively. However, having identified an integrable long-range spin 
chain which properly describes the dilatation operator up to three loops, we should take 
it seriously. Integrability is a subtle and tight structure, and not found abundantly. 

We have not been able to explicitly mention in this introduction many interesting, 
closely related recent contributions to this subject of study, see in particular [37] . 

2 The Inozemtsev spin chain and J\f = 4 gauge theory 

2.1 Long-range spin chains with elliptic and hyperbolic inter- 
actions 



More than one decade ago, Inozemtsev j22j proposed a spin chain with long range inter- 
action which interpolates between the Heisenberg and Haldane-Shastry c 32- spin chains. 
His Hamiltonian takes the form 

L L-l 

^EE^aW^-^)' ( 2 ) 

j=l n=l 

where the operator Pj k permutes the spins at sites i and j, and the interaction strength 
is defined in terms of the elliptic Weierstrass function "P£ j7r / re (n) with periods L and in j 'k. 
We recall that the Weierstrass function can be defined by the following series 

L,7r / K z 2 ^-^ \(z — mL — inir / k) 2 (mL + inir / k) 2 J ' 

where the prime means that the term with m = n = is omitted. When one of the two 
periods becomes infinite, the Weierstrass function becomes a trigonometric or hyperbolic 
function. If k — > 0, this function goes smoothly to the Haldane-Shastry interaction, 



HmV L /K (z) = [j i i — - 



while for L — > oo the interaction takes the form of a hyperbolic function, and decays 
exponentially with the distance, 



lim Vl,-k/k{z) = ( — g — + - j 
l^oo \ smh kz 6/ 



The Heisenberg limit can be recovered in the limit k —>■ oo, where the imaginary 
period vanishes. To be able to properly take this limit we define the "coupling constant" 



and we use a redefined interaction strength, by removing an additive and a multiplicative 
constant 



The Heisenberg interaction can be obtained to leading order in t of the Inozemtsev 
interaction 



Going beyond the first order in t, we notice that ht^n) starts as t n (plus higher orders 
in t), as long as n < L/2, therefore the Inozemtsev interaction involves spins separated 
by a distance of at most n at order n in t. This property agrees with the structure of 
the perturbative dilatation operator in the planar limit |24j . 

Inozemtsev performed a detailed study of the Hamiltonian (J2J) (see [33] for a re- 
view) and gave convincing arguments in the favor of its integrability, by finding the 
corresponding Lax pair 5 . In addition, he explicitly constructed the first of the higher 
charges commuting with the Hamiltonian. The existence of this first higher charge usu- 
ally indicates complete integrability for spin chains. However, the general form of the 
commuting charges or the monodromy matrix are still out of reach. The connection 
with the Calogero-Moser Hamiltonian with elliptic potential presented in [33] may help 
finding the exact solution of this chain. Among the aforementioned limiting cases, the 
Heisenberg case k — > oo and the Haldane-Shastry — ► are completely solved, the 
first using the Bethe ansatz and the second by exploiting the Yangian symmetry [3§] . 

In the following, we are going to concentrate on the hyperbolic case, where the length 
of the chain becomes infinite (i.e. L — > oo), and the interaction takes the simpler form 



that is, the periodicity effects of the interaction are negligible. This case is easier to 
treat analytically, but still quite interesting from the point of view of comparison with 
the known results for the M = 4 SYM theory. In particular, it allows to study BMN 
operators at finite (but large) J ~ L. The full elliptic potential takes into account the 
effects related to the periodicity of the chain, therefore one could hope that it contains 
the information concerning "short" operators such as Konishi. 

In the hyperbolic regime we can use an analytic continuation for the monodromy 
matrix and for the conserved charges of the Haldane-Shastry spin chain [313 IHH l4*Tj . 
The Haldane-Shastry conserved charges still commute (and commute with the Yangian 

5 According to Inozemtsev jjH], who invokes a result by Krichever |38j . this is the most general 
Lax pair corresponding to a spin chain with long range interaction and pairwise exchange given by 
permutation operator. However, his arguments do not directly apply to the case when the interaction 
involves more than two spins, and we should keep an open mind as far as concerns the existence of 
integrable long range spin chains of a more general type than Inozemtsev. For these putative chains, 
the subtraction of higher charges from the lowest charge (i.e. the Hamiltonian) would not allow one to 
eliminate interactions involving more than two spins. 




m=— oo 
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generators) if one replaces the coordinates Zf- = e 2mk ^ L by Zk = t k [10] 6 - We recall that 
the Yangian generators are given by 



3 

no „-_a6c \^ Z j j~ Z k J, ^ c t A \ 



Zj 



Let us emphasize that the Yangian is a true symmetry (i.e. it commutes with the con- 
served charges) only for the hyperbolic chain, when the length of the chain is infinite. 
Otherwise, like in the case of Heisenberg spin chain jl2], in the full elliptic case the 
Yangian symmetry is broken by "boundary terms" . The higher conserved charges of the 
hyperbolic spin model are also obtained by simply replacing Zk = t~ k in the expression 
ofor the Haldane-Shastry conserved charges [3§1 HT] 



E 



Zj ^-(Pij - 1) , 



ZijZji 
ij J J 



. . , ZijZjkZki 
IJK J J 

h = J2 f iZjZkZl - i) - 2 E ( 7^-) 2 - l ) » ( 5 ) 

■1.1 Z ij Z jk z kl^li ■■ \ z ij z ji J 

where the prime indicates that the sum is over distinct summation indices, Zij = Zi — Zj, 
and Pijk, etc. represent cyclic permutations. A systematic procedure to find all the 
conserved charges was given in [H], although the explicit construction of these charges 
becomes more and more involved at higher orders. 

The diagonalization of the conserved charges is achieved, in the case of nearest- 
neighbor interaction, by the algebraic Bethe ansatz. In the case of the long-range inter- 
action, the algebraic Bethe ansatz does not work any more, just like in the case of the 
Haldane-Shastry spin chain. There, the diagonalization can be performed by using the 
Yangian symmetry and the relation to the Calogero- Sutherland model [SHIES!- In the 
hyperbolic case it is still possible to use the so-called asymptotic Bethe ansatz, valid in 
the regime where L is very large but still finite. The phase acquired by a quasi-particle 
(magnon) traveling around the circumference of the chain is given by the phase gained 
by its scattering with all the other magnons in the system 



exp (ipjL) = exp 



( M \ 



(6) 



6 Thc generating function for the conserved charges of the Haldane-Shastry Hamiltonian was identified 
in |41| by taking the limit A — > oo of the monodromy matrix for the Calogero-Sutherland model with 
spin. The same argument can be easily applied to the hyperbolic case, since in this case the leading 
term of the quantum determinant is trivial, Ao(u) = 0, as well as its derivatives. It follows that the 
next-to-leading term Ai(u) in A -1 in the expansion of the quantum determinant can be used to generate 
the commuting charges in exactly the same fashion as for Haldane-Shastry spin chain. 
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where xiPjiPk) is the phase shift for the scattering of two magnons of momenta pj and 
Pk and M is the number of magnons. Inozemtsev computed the phase shift for the 
scattering of two magnons, 

cot = tpipj) - <p(p k ) , 

^ = 2^ Ci (I) - i Ci (!) ' (7) 

as well as the magnon energy, 

S(P) = ^C. ( f\ - M ( f] - S(P) , (8) 
m k \2k j 4«r \2k j 

where V\{z) = d(i(z)/dz denotes the Weierstrass elliptic functions of periods 1 and itt/k, 
and the quasi-periodic zeta function is odd, Ci( — z ) = ~Ci( z )- The function ip{p) is also 
quasi-periodic, ip{p + 2n) = <p(p) and ip(p + 2in) = Lp(p) — i. While the phase shift and 
magnon energy ()7|8)1 are exact, the Bethe ansatz formula (jH} is asymptotic, i.e. it is valid 
for chains of finite but large length L. It is easy to check that in the limit k — > oo, the 
magnon energy and the phase shift reduce to the corresponding Heisenberg values, in 
particular <p(p) — > (1/2) cot(p/2). 



2.2 Comparing the conserved charges of the Inozemtsev spin 
chain to the dilatation operator 

If we assume that the dilatation operator in the su(2) sub-sector of M = 4 theory is 
integrable to all orders in perturbation theory, as first conjectured in it should 
correspond to some integrable long-range spin chain. We believe that the Inozemtsev 
spin chain is an interesting and natural candidate. As we shall prove, it is definitely 
capable of describing the dilatation operator in the subsector to at least three loops. 
However, if the describtion extends to four loops and beyond, perturbative BMN scaling 
would have to be violated at the four loop level. 

If the dilatation operator contains, beyond two loops, interaction terms which involve 
more than two spins, the Inozemtsev Hamiltonian is not, by itself, sufficient to reproduce 
the former 7 . In this case, one should consider a linear combination of the even parity 
conserved charges of the Inozemtsev spin chain 8 . The most general expression for the 
dilatation operator will then be 

D(X)=L + f 1 (X) I 2 (t) + f 2 (X)h(t) + f 3 (X)h(t) + ... (9) 

7 The three-loop dilatation operator proposed in |27j does contain four-spin terms, as these are nec- 
essary for perturbative BMN scaling. If we drop this requirement, the yet undetermined two constants 
29 discussed in the introduction may be chosen so as to eliminate the four spin terms. The future will 
tell. 

8 In |27j it was argued that only connected planar diagrams contribute to the dilatation operator 
(see also A completely rigorous all-orders proof is not easy due to the known problems with 

regularization by dimensional reduction at large loop order in Af — 4 gauge theory. Non-connected 
diagrams would require non-linear combinations of the conserved charges. 
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where t = t(X) is such that \im\^ t/\ = 1/167T 2 , /fe(A) are functions regular at A = 
and l2k{t) are the even parity conserved charges of Inozemtsev spin chain. 

To determine the precise relation between the 't Hooft coupling constant A and the 
parameter t, as well as the functions /fc(A), we have to use arguments going beyond 
perturbation theory. The first example of such an argument concerns the behavior of 
the rapidity (p(p) in the BMN limit, where the coupling constant A and the chain length 
L are both large, A = A' L 2 . Since it enters the Bethe ansatz equations, and therefore all 
the observables, the combination (p(p)/L should have a well-behaved BMN limit. The 
rapidity (JJj) can be written as an infinite sum 



1 

- cot 
2 2 



V ly 

9 9 ^ 



, cot(- — inn) + cot(— + inn) 
2^L y 2 1 y 2 ' 

n>0 



(10) 



- cot - + 2 V — 
2 2 ^ 1 

n>0 v 



t n sin p 



t n ) 2 + 4t n sm 2 (p/2) 



where in the second line we have used the perturbative parameter t instead of the pa- 
rameter k. Expanding this result to the first few orders in p we obtain 



cp(p) = ^cot| + 2pJ2jl 



n>0 



t 



n\2 



3 ^ 



n>0 



t n (l + 4t n + t 

(i - t n y 



2n\ 



+ 



In the BMN limit, the leading behavior of the magnon momentum is p ~ and 
therefore the rapidity scales as tp ~ 1/p ~ L. In order to obtain a finite but non-zero 
result for the rescaled rapidity tp/L in the perturbative regime A ~ where k — ■> oo 
(i.e. t — > 0) we see that we need to identify 



A 



16tt 2 



n>0 v 



- t 



n\2 



E 

n>0 



1 



4sinh (kti) 



On then finds that 



L 



1 2 A 



3L 4 



167T^ 



+ 6 



A 



+ 



A' 



8 



7T- 



(Lp) 



k'2 



1287T' 



-{Lpf + L 1 



16vr 2 
A' 3 



12 



1024vr 6 



167T 2 



(Lp) 3 + . . . 



as 



(Lp) 3 + ... 



oo (12) 



where we have used the finite BMN coupling constant A' = X/L 2 . 

We see that at C(A' 3 ) a quadratic divergence appears in the perturbative expression 
for the rescaled rapidity f/L. The consequences of this perturbative scaling violation 
will be discussed in section 2.4 below. 

It is interesting to see that a finite rescaled rapidity ip/L may be obtained provided we 
stay in the strong-coupling regime where k — > (i.e. t — *■ 1) and thus A = 2tt a /3k 2 — > oo. 
Here, it is useful to use a dual representation for the rapidity ip(p) 



VIP) = — coth 1 > 

YyF > 2k 2k 2k ^ 



m>0 



coth ; ^ + coth - u 



2k 



2k 
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where the last term insures the desired quasi-periodicity properties, <p(p + 2ix) = ip(p) 
and (f(p + Hk) = ip(p) — i. The terms with m > form a power series in the dual 
parameter e~ 2n l K and they vanish exponentially in the BMN limit, since in this case 
1/k ~ L — > oo. The remaining part obviously obeys BMN scaling and may be written 

as 

v(p) = l\f¥ ,co Kj¥'' (14 > 

where we have replaced <p(p)/L by <p(p) and pL by p. 

In principle, the explicit knowledge of all the conserved charges of the Inozemtsev 
chain and of their eigenvalues would allow to obtain information about the functions 
/fc(A) in equation ©. Since, at least for the moment, we lack this knowledge we have 
to restrict ourselves to comparing the relation (jHJ) to the perturbative results which 
are already known t 27 , [29J. Let us then compare the expression (JOJ) to the dilatation 
operator, 



fc>0 

as stated in [27]. Rewritten using commuting permutations, the contributions to the 
dilatation operator up to three loops are 

D = L, (15) 
D 2 = 2L-2^P M+1 , 

i 

D 4 = _ 6jL + 8 ^p.. +1 _2^P M+2) 

i i 

Dq = 40L — 56 Pi,i+i + 16 Pi,i+2 — 4 (-Pj,i+3-Pi+i,j+2 — ^,4+2-^+1,1+3) • 

i i i 

We recall that D 2 was obtained by diagrammatic computation in M = 4 SYM theory 
and the expression of D4 was checked in several ways |21],j2Sl- The expressions for 
D$ were obtained initially |27| HI] supposing perturbative integrability and imposing 
perturbative BMN scaling. The expression for D$ was subject to controversies, since it 
disagrees with the results from the near plane- wave string theory [2H] • Later on, Beisert 
checked that the integrability of D Q is insured by supersymmetry, and that the two 
free coefficients are fixed if one imposes perturbative BMN scaling. 

The conserved charges of the Inozemtsev spin chain can be expanded in series of t, 
and therefore of the coupling constant A. By inverting the series (jllj) we obtain 

A / A V . / A 



i= 16^- 3 U^ +1 Hl6^J + "" (16) 
Inserting this expression into (J3J and expanding in powers of A we get 



n>0 



( 76^ ) ^ 
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The explicit expression for the first few orders in the Hamiltonian is 

J 2 ^ = 2L — 2 ^ ^ Pi,i+l j 
i 

^2 = ^ P«,i+1 — 2 Pj,i+2 , 
j i 

i^P = ~iq -Pj.j+i + 12 -Pj,i+2 ~ 2 -Pj.i+3 ■ 

i i j 

In the following, we use a redefinition of the fourth conserved charge, /4(A) = /4(A) — 
g^a /2(A). This is done for the purpose of having the fourth conserved charge of the 
Heisenberg chain as the first non-trivial order in the A expansion 

j(l) _ j(2) _ q 
M — J 4 ~~ u > 

if = -8L + 8^P M+1 -4^P M+2 + 4^P M+3 

i i i 

— 8 (-Pj,i+3-Pi+l,i+2 ~~ -Pj,i+2-Pi+l,i+3) 5 

i 

It is not obvious that a similar redefinition can be done for all the conserved charges, 
but at least we know that we can arrange to have lg = 0. 

We can now reconstitute the perturbative dilatation operator (fTKJ) using the Inozemt- 
sev conserved charges 

D, = I«, 

At = 4 2) -3/«, 

£ 6 = 4 3) + - 3/f + 22/^ . (17) 

If we reinstate the constants which were fixed by imposing perturbative BMN scaling 
j2Tj , the most general form for the third loop dilatation operator allowed by perturbative 
integrability is 

D & = 4 3) + ajf ] - 3/< 2) + a 2 J« . (18) 

Comparing to the expression JOJ), we conclude that /i(A) = 1 — (2ai + 3)y^2 + • • • and 
/2(A) = CL\ + . . ., the constant a 2 appearing at next order in A. In conclusion, we can 
reproduce the dilatation operator up to three loops using the conserved charges of the 
Inozemtsev spin chain. 

2.3 Perturbative results from the asymptotic Bethe ansatz. Di- 
mensions of the BMN operators 

In the following, we are going to compute perturbatively the dimensions of the BMN 
operators up to third order, supposing that the dilatation operator is integrable and its 
integrability is governed by the Inozemtsev spin chain. The strategy is to compute the 
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magnon momenta by using the asymptotic Bethe ansatz equation (JBJ) and to substitute 
them into the expression of the dilatation operator known from perturbative calculations. 
The asymptotic Bethe ansatz equations have the form which we can rewrite as 

e** L = f] <£i^£h±i , (19) 

where M is the magnon number and the rapidity (fj = <p(pj) is given by the expression 
(jlOj) . To see the connection with the usual Bethe ansatz for the Heisenberg chain, we 
can replace the momentum by a rapidity variable Uj = | cot y such that the left hand 
side of eq. llQj) becomes 

e i PjL = ( Ui + i/V ' 
\u j -i/2 / 

while it is clear from eq. ljlOj) that (fj = Uj + 0(t). More precisely, up to third order in t 
we have 

<p(jp) = u(p) + 2t sinp + 2t 2 sinp (1 + 2 cosp) + 8t 3 sinpcos 2 p + . . . . 

We can perturbatively solve the Bethe ansatz equation for the two-magnon case, with 
the two magnon momenta satisfying p± = —p2 = p. Since we will compare our results to 
the anomalous dimensions of the BMN operators 

TrZ J $ 2 + ... , 

we are going to use J + 2 to denote the length of the chain, 

L = J +2 . 

The Heisenberg solution for the magnon momentum is 

2-7m 

P = Pn = 



J+l ' 

and the corrections to this value are given by the equation 

8tsinpsin 2 f 8t 2 sin p sin 2 f / .p 2 p 
P = Pn — x 7TT^( 8cos 2" 4C ° S 2 + 1 

32t 3 sinp sin 2 f / p 6 p 4 p 2 p 

- ' 64 cos - — 96 cos - + 48 cos - — 16 cos - + 3 + . . . 



3(J+1) V 2 2 2 2 
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The perturbative solution of the Bethe ansatz equation up to third order in t is 
St sin p n sin 2 ^ 



P = Pn' 
+ 



+ 



J+l 

64t 2 sinp n sin 4 ^ (4 cos 2 ^ - l) 8t 2 sinp„ sin 2 2& (8 cos 4 2& - 4 cos 2 £& - l) 

(J+l) 2 7+1 
512t 3 sinp n sin 6 (24 cos 4 ^ - 13 cos 2 ^ + l) 

128 1 3 sin p n sin 4 (^) (48 cos 6 ^ - 44 cos 4 ^ + 4 cos 2 2* + l) 

(J+lJ 2 

32t 3 sinp n sin 2 2* (64 cos 8 ^ - 96 cos 6 ^ + 48 cos 4 ^ - 16 cos 2 *p + 3) 

3(7+1) 



+ 



Using the perturbative expansion (|16|) for £ in powers of A, we can transform this expres- 
sion into a series in A. Then, we insert the solution for the momentum in the eigenvalue 
formula for the dilatation operator, which is a consequence of (|15)1. 

£(a) = j + 2 (i + A sin2 p _ sin 4 £ + sin6 P 

v ; V 2vr 2 2 8tt 4 2 16vr 6 2 

Using Maple to perform the resulting series expansion, we checked that the two-magnon 
energy at order k in A, E^ k \ has the behavior predicted in [27], up to order k = 3 

= ^sin 2 ^, (20) 

77.(2) 1 ■ 4 Pn / 1 , 1 2 Pn 

= -^ sm tU + 7+T cos T 

1 Pn /l I cos 4 g + f cos 2 I § cos 4 f - I cos 2 g 

E = ^ sm YU 7+1 + (TTIp 

The agreement with the perturbative results in gauge theory is perfect, up to three loops, 
as far as the anomalous dimensions of the large J BMN operators are concerned. For 
completeness, let us look at one more order and, in particular, to compare with the 
conjectured results of |Hj. 



2.4 Fourth order dilatation operator and the breakdown of per- 
turbative BMN scaling 

The four loop dilatation of M = 4 gauge theory is not yet rigorously known. In [44 
a proposal for its structure was put forward, based (1) on extending the planar pair 
symmetry, due to integrability |2?J, and (2) on insisting on perturbative BMN scaling. 
As we shall see below, this proposal will turn out to be inconsistent with the Inozemtsev 
integrable structure. We suspect that the constraints coming from all-order integrability 
become tighter at higher orders and that it will become impossible to satisfy at once 
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integrability and perturbative BMN scaling. However, this does not mean automatically 
that [33] is wrong. Indeed, it is not excluded that several integrable long-range spin 
chains exist (see the comments in footnotes 4 and 5). 
The proposal in 03] can be rewritten as 



D 8 = 222L- 648^P M+1 + ^^P M+2 - ™^P M+3 + 6^P M+4 (21) 

i i i i 

+ P,i+lP+2,j+3 H g — P,i+3p+l,i+2 — 366 Pj,t+2p+l,i+3 

i i i 

+ ~ y^(Pi,t+3p+2,i+4 + Pi,i+2pi+l,j+4 — Pj,i+4Pi+2,i+3 — Pi,i+4pi+l,i+2 
i 

~ Pi,i+4:Pi+l,i+3 + P,i+3p+l,i+4) 5 

where we have chosen the value j3 = —428/3 in the expression of [33], since this value 
allows to reconstitute the term in parenthesis, which is a very natural combination of 
permutation. We recall that (3 corresponds in [44J to a unitary transformation of D 8 and 
that it does not affect the eigenvalue of D 8 . 

We also expand the first two even Inozemtsev conserved charges to fourth order in A, 

1*2 — 62 ^ Pi,i+i — 78 ^ Pj,i+2 + 18 ^ P,i+3 — 2 ^ Pj,j+4 , 

i i i i 

5 4) = 48L-36^P M+1 + 16^P M+2 -40^P, l+3 + 12^P ii+4 

i i i i 

- 8^P M+1 P i+2)i+3 + 56 Pi,i+3p+l,j+2 — 48 Pi,i+2Pj+l,i+3 
i i i 

+ 8 {Pj,i+3Pi+2,i+A + Pi,i+2pi+l,i+4 — Pi,j+4pi+2,i+3 ~ Pi,i+4p;+l,i+2 
% 

~ Pi,i+4pi+l,i+3 + P,i+3p+l,i+4) • 

We see that ansatz ()21)1 can be written as 

2 ~ 

D 8 = I 2 + q-^4 4 '' + lower order terms , (22) 

where "lower order terms" stands for a term of the form ^j(P,j+iP+ 2 ,i+3 + P ) i+3p+i,i+ 2 ) 
which we expect to come from plus a linear combination of J 2 , I 2 , J 2 ■ 

What we learn from this exercise is that, imposing the integrability perturbatively, 
as it was done in |2"7] IHj . we obtain the building blocks of the Inozemtsev spin chain 
(modulo the contribution of the sixth conserved charge, for which we do not have a 
quantitative estimate). However, the coefficients multiplying in (|T7|) disagree with 
the all-order expression (jH]). The expression (jH]) implies that the coefficient of 1^} in P 2n 
should not depend on n. This is not the case for 1^ and 1^ appearing in P 6 and D 8 in 
equations (fT7|). (jUJ). In order to be compatible with the proposal (JHJ), the fourth order 
dilatation operator has to be 

D 8 = I? - 3if + oaif + 6x J« + a^f) + 6, jf + bjf . (23) 
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A possible explanation of the mismatch is related to the fact that the coefficients in D§ 
and D s which were left free by perturbative integrability were fixed imposing perturbative 
BMN scaling |27| 144 j. However, in (|23|) there are three free constants left, once we have 
fixed, at third order, a\ and a 2 , while Beisert obtains |S] six constants to be fixed by 
perturbative BMN. This point certainly deserves a better understanding. 

Independently of the previous analysis, we know that the Inozemtsev spin chain 
cannot reproduce, starting from the four-loop order, an anomalous dimension which 
obeys BMN scaling. The reason is our finding in eq. (jl2|) : The rescaled rapidity (p(p)/L 
does not have an expansion in A consistent with perturbative BMN scaling. The first 
violation appears at third order in A, and therefore will contribute to the correction to 
the two magnon energy at the fourth order. 

Although the rapidity (fTT)|) does not obey perturbative BMN scaling (when A is small), 
it obeys BMN scaling in the BMN regime, when A is large. Technically, we see that the 
perturbative BMN scaling is spoiled by terms of the type e~^, which vanish in the BMN 
regime but have a well-defined, non-trivial expansion around A = 0. 

3 Stringing spins and spinning strings at two and 
three loops 

3.1 Generalities 

In the last chapter we showed that the Bethe ansatz for the Inozemtsev chain is indeed 
capable of reproducing the correct two- and three-loop anomalous dimensions of finite 
length two-impurity operators. Let us now turn to the "thermodynamic" situation of 
"long" operators with many impurities; i.e. we will study the operators of "length" 
L = Ji + J 2 and "magnon-number" M = J 2 

TiZ Jl $ j2 + ... , (24) 

in the limit where both J\ and J 2 are large. At the one-loop level anomalous dimensions 
for these operators were obtained in [IB] and it was established in [T%| HH] that the result 
agreed with string theory predictions of Arutyunov, Frolov, Russo and Tseytlin Hl)| 
ITTj . in two cases: The folded string (corresponding to the ground state) and the circular 
string (corresponding to an excited state). As we explained above, the Inozemtsev chain 
turns out to be intrinsically inconsistent with perturbative BMN scaling, starting at 
four loops. Just like BMN, the all-loop proposal of Frolov and Tseytlin also predicts an 
expansion in the BMN coupling A' = A/ J 2 . We immediately infer that the perturbative 
Inozemtsev ansatz therefore has no chance of reproducing these predictions beyond the 
three-loop level. On the other hand, up to three loops BMN scaling is possible, and, at the 
two-loop level, proven. As was shown in the previous chapters, the three-loop dilatation 
operator of [27J can be emulated by the Inozemtsev chain. It is therefore clearly very 
interesting to work out the Inozemtsev predictions at two and three loops and compare 
to string theory. This will be done in the next sections, for both the folded and the 
circular string. At the two-loop level we can be sure that we are indeed computing the 
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perturbative gauge theory result. This is in contradistinction to the three-loop case, 
where we need to assume the absence of the terms that break perturbative BMN scaling. 
Let us recall the string result for the two relevant cases of spinning folded and circular 
strings. It may be written in the following concise way. For the folded string we have 
the parametric expression 



J2 V f Ji V 4 



(25) 



K(t) -E(t) J \E(t) J 7T 2 ' 

t , (26) 



£ V ( Ji V 4 



K(t) J \E(t) J it 2 
while the circular string energy £ is given by 

Ji V ( J2 V 4 



K(t) - E(t) J t)K(t) - E(t) J ttH 



(27) 



where £,J\,J2 are the string energies E and the angular momenta J\,Ji in units of 
the effective string tension A (alias square root of the 't Hooft coupling), E = y/\ 8, 
J\,2 = vA Ji,2- K(t) and E(t) are elliptic integrals of, respectively, the first and second 
kind. In both cases the upper equations (|23|) . (|27j) determine the parameter t 9 as a 
function of J\,J2, while E is then obtained by substitution of t into, respectively, the 
lower equations (|25 jl .(pS jl . These expressions may be expanded for large J = J\ + J2] 
this results in an expansion of the energy as a Taylor series in the BMN coupling A' 
through the identification (cf. eq.((T}) 

^ = j~ 2 - (29) 

In order to obtain a finite result for the energy £ = J + . . . we should introduce a rescaled 
string energy (=lowest charge) Q2. 

Q2 = j. (30) 

One then has 

00 

Q 2 = 1 +E Q{ " ] ( A ') fe ' ( 31 ) 

k=l 

where the expansion coefficients only depend on the dimensionless ratio (termed 
"filling fraction" below) 

— ^2 J2 , qo n 

a =j=r (32) 

9 The parameter t is not to be confused with the Inozemtsev coupling constant t = e~ 2K of the 
previous chapter. 
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with J = J i + J 2- Since the expansion eq. fjUlj) is in integer powers of A', the basic pro- 

(k) 

posal of Frolov and Tseytlin has been that the k-ih expansion coefficient Q 2 should 
be reproducible by a fc-loop perturbative gauge theory calculation of the anomalous di- 
mension of the operators eq. (|24|) . This proposal generalizes the one of BMN, as the case 
of "J2-impurity= J 2 -magnon" BMN operators may be recovered from the folded string 
solution eqs.(J2HJ),(j2H|) in the limit where J 2 stays finite (i.e. J 2 — > 0), but the length of 
the operator still goes to infinity: J ~ J\ — > 00. In this case eq.(jSH) becomes 

Q 2 = 1 + y (v/TTA 7 - l) (33) 

(k) 

i. e. the Q 2 become proportional to the Catalan numbers: 

Q * -7 8( "4 } k\(k-i)\' (34) 

As discussed in the introduction, both the BMN and the Frolov-Tseytlin proposal make 
the highly non-trivial assumption that in the AdS/CFT system two different scaling 
procedures are nevertheless equivalent: Recall that in perturbative gauge theory the 
't Hooft coupling A is small, while the string calculation requires A to be large. We are 
now ready to test the Frolov-Tseytlin proposal beyond the one-loop level, with the help 
of our long range Bethe ansatz eq.([T§ J) ,, ([Tu ft - 
s' Thermodynamic limit of the Inozemtsev-Bethe equations 

The thermodynamic limit of the Inozemtsev-Bethe equations is obtained in close analogy 
to the procedure employed for the Heisenberg model, cf. [IB]- One takes the logarithm 
of eq.(Hnj) and assumes the length L = J = J\ + J 2 to be large. This gives 

J 2 1 

J Pi = 2vm 7 - + 2 V . (35) 

In the case of the one- loop Bethe ansatz one works with the Bethe roots Uj, or, equiva- 
lently, the momenta pj. In the present case one has, in addition to the {pj}, a second set 
of discrete variables, namely the rapidities {<£>•,}. The two sets are related through the 
coupling constant dependent relation eq.(fTU|). In order to obtain equations which stay 
close in form to the original Bethe ansatz it is most convenient to use eqs. ()35|) to deter- 
mine the distribution of the rapidities {y?j}, instead of working with the momenta {pj}- 
We can think of the rapidities as "deformed" Bethe roots, since we have ipj = Uj + 0(A). 

In the limit where the magnon number M = J 2 gets large, the rapidities are expected, 
after a rescaling, to accumulate along smooth contours in the complex (p plane, just as 
its undeformed leading order values Uj 

^-v with ^) = 7E 5 (^t)> ( 36 ) 
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where p(ip) is a distribution density which is normalized to the filling fraction a = J2/J, 



Jc 

and C is the support of the density, i.e. the union of contours along which the rapidities 
are distributed. The Bethe equations (|33J) then turn into singular integral equations: 



nc(<p) is the mode number at point p. It should be constant along each contour. Here and 
in the following the slash through the integral sign implies a principal part prescription. 
This equation generalizes the continuum Bethe equations of [18J. On the right side the 
non-trivial potential p p(ip)/2 appears, which is given by the scaling limit of the inverse 
of the rapidity p(p) p/2 in eq. ljlO)) . Eq. (|38)l should be complemented by the momentum 
conservation condition j^, dip p(<p) p(<p) = and the constraint dtp p((p) nc( v ) = 0. 

We would like to stress that the Inozemtsev-Bethe equation eq.(j2HJ) is expected to 
determine the spectrum of the long-range spin chain (and therefore, potentially, of the 
operators in eq. (|24|) ) quite generally in the regime J±, J 2 large (neglecting, however, the 
"wrapping" modes, as explained before). It is not necessarily restricted to perturbation 
theory. The information on the observables (namely, the commuting charges) of the long- 
range spin chain should then be contained in the moments Qk of the rapidity distribution, 
in generalization of the result of |21j.|2Uj: 



However, as far as determining the energy (=lowest charge) of the spin chain is 
concerned, finding the distribution of rapidities is not yet everything. We also need 
to uncover the relation between the "natural" observables Qk and the observables we 
would like to measure, i.e. in particular, the eigenvalue of the dilatation operator. In the 
perturbative regime, up to three loops, this problem was solved in the last chapter. To 
this order, we are now ready to obtain the anomalous dimension of the operators eq. (|2"4*|) 
in the two cases of interest (folded and circular). 

3.3 Folded case 

The folded case corresponds to the ground state of the operators eq.(J2H). The qualitative 
rapidity distribution is expected to be as in the one- loop case [TH]. We have precisely 
two cuts symmetrically placed to both sides of the imaginary axis, with, respectively, 
mode numbers ±n. Using the perturbative result of the last chapter, we find (putting 
the mode number n to one), up to three loops, the integral equation 




(37) 




(38) 




(39) 




<p' 2 - <p 2 4 2 




A' A' 2 



(40) 
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which may be treated in a fashion identical to the one appearing in the multicritical 
O(n) matrix model [IS]. The constants ci,c 2 are found from the results of the previous 
chapter to be 

I 3 

and C2 = — r- (41) 



327T 2 5127T 4 

Let us define the moments Q 2 k of the ip distribution by 



Q 



2k 



y (p2k 



(42) 



The energy Q 2 {}■£■ the anomalous dimensions divided by J) is then determined 10 by 

A' 



Q'2 



47T 5 



Q 2 - e x A' 2 Q 4 - e 2 A' 3 Q 6 , 



where the constants ex, e 2 were found in the last chapter 

3 , 5 



ei 



64tH 



and 



('2 



5127T 6 ' 



(43) 



(44) 



It is easily verified that these equations are a generalization of the one-loop ones presented 
in ^Hj; their solution proceeds as before and we will thus be rather brief; for further 
details, see [21]. The normalization of the density is 



dip p{ip) 



a 



(45) 



with the by now standard trick [T%] IT9"] |2*T] of analytic continuation to (formally) negative 
J, resulting in a negative filling fraction a = J 2 / J < (this flips the complex cuts onto 
the real axis). It is again useful to introduce a generating function for the moments Q 2 k- 



hk 



2k 



i.e. 



k=l 



dip' p(tp') 



12 



If' 2 — (p 2 



(46) 



We find the following generalization of the result eq.(2.24) in |21j : 



H{p) 



a 
2 



A' 

+ T + Ci —r + C 2 



1 

4 



'b 2 



f 9 



a 2 — ip 2 



n 



ip 2 )(a 2 — ip 2 



ci A' — + c 2 A 

cp z 



12 



-Q 



~2 2' 9 

a A — (p z 



a 2 + b 2 1 
2a 2 6 2 <^ 2 



(47) 



10 Formally the moments Q2fc, and in consequence also the energy Q2, have an infinite regular expansion 
in the BMN coupling A'. The leading order 0(X'°) contribution to Q^k is nothing but the 2fc'th one-loop 
commuting charge, cf. Obviously we should discard all terms higher then C(A' 3 ) in Q2 since the 

thermodynamic Inozemtsev-Bethe equation ea. H4()(l is accurately describing perturbative gauge theory 
only to at most C(A' 3 ). As we discussed, in light of the results of (2Hj, we can actually only be completely 
sure about this validity to C(A' 2 ) and need to assume the vanishing of two so far undetermined constants 
to ensure BMN scaling and, in consequence, validity of eq.lQUJ) to C(A' 3 ). 
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where the modulus q is, as before, q = 1 — a 2 /b 2 . The II function is the elliptic integral 
of the third kind; for our conventions we refer to [21]. One then finds that the left edge 
a of the cut is determined through the equation 



\-K(q)a + c 1 A' (1-f)^ 

while the right edge b is, of course, b = a/y/T 
the filling fraction to the modulus q: 

1 2aE(g) 

2 " 



+ c 2 A' 2 (1 - q 



0. 



(48) 



q. Finally, one has an equation relating 



a = - 



c\ A 



+ c 2 A 



12 



(2-q) x /T^ 



(49) 



K(q) and E(g) in eqs. (f4*Kj) . (|4T)|) are elliptic integrals of the, respectively, first and sec- 
ond kind. We are now able to obtain the gauge theory result up to three loops, using 
eq.(jHH).This is done by expanding both the modulus q = q Q + q\\' + q 2 \' 2 and the in- 
terval boundary a = a$ + a\\' + a 2 A' 2 in powers of A', and eliminating with the help 
of eqs. (|4*H]) . f|4*n|) all constants except for the leading order modulus go, which, as before, 
conveniently parametrizes the filling fraction a through 

h 1 E( 9o ) 1 

This gives the final, three-loop result 

1 



2VT^ K(q 



(50) 



I + _K(g )(2E(g ) 



(2 - go)K(go)) A' 



(51) 



1 



8vr 4 
1 



K(g ) 3 (4(2 



go)E(go) - 



8q + q 2 )K(q )) A 



12 



4tt 6 (E(<zo) - K( go ))(E(g ) - (1 - ?o)K(g )) 



8g + 3gg)E(g 



(2 - q )(12 - 12g + q 2 )E(q ) 2 K(q ) 



+3(1 



q )(8 - 8q + q 2 )E(qo)K(q ) 2 
-4(l- qo ) 2 (2-q )K( qo ) 3 



A 



/.3 



We notice that the two-loop (0(X' 2 ) contribution in eq. (|51jl very closely resembles the 
eigenvalue of the fourth higher charge Q±, cf. eq.(2.25) in |21j . This finds a natural 
explanation in terms of the structure of the energy expression eq. pH|) : Indeed the two- 
loop energy is "mostly" given by the contribution from the fourth conserved charge of 
the Heisenberg XXX chain, with a "small" correction, namely the A' correction to the 
second Heisenberg charge. 

Let us now compare this to the string prediction eqs.()25 )) .([26 p . The latter is conve- 
niently rewritten in Lagrange-inversion form, using eqs. ()29|) . (j32j) : 

2 / T71 / _L \ \ 2" 



A' 



Q 2 



71 



4K(t 



K(t ) - E(t ; 

K(t) - E(t) 



V E(t) J 



(52) 




t) 



f E(to) 



Kgo) - ggoj 
K(t) - E(t) 
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where the "string modulus" to is found to parametrize the filling fraction through 

a = 1 ~WhY (53) 

One now expands both equations in t around the point to- The upper equation gives A' 
as a power series in the variable (t — to). Power series inversion yields t as a power series 
in A'. Finally, substitution of this series into the expanded lower equation gives Q2 as a 
series in A', i.e. the sought coefficients in eq. lpTT]) . We then find, up to third order in A', 

Q 2 = 1 + 4cK(i ) [E(t ) - (1 - io)K(to)] A' (54) 
+^K(t ) 3 [(1 - 2to)E(to) - (1 - t ) 2 K(t )] A' 2 

7T 

4 K(t f 

+ vr 6 E(t ) 2 - 2(1 - to)E(to)K(i ) + (1 - to)K(to) 2 X 

x ((1 - 7t + 7t 2 )E(t f - (1 - t )(3 - Ut + 7t 2 )E(t ) 2 K(to) 

+(1 - t )(3 - 2t )(l - 3t + t 2 )E(t )K(t ) 2 - (1 - t ) 4 K(t ) 3 ) A' 3 

+0(A' 4 ), 

In order to compare to gauge theory, we need to parametrize both results in the same 
way, cf. eqs. (l5U]) . (|53Jl . In fact, we can apply the Gauss-Landen transformation of jTHj, 
which relates the string and gauge theory modulus through 



t0 = Jl^^L. (55) 

to rewrite the three- loop string energy Q2 in terms of the gauge modulus qo'. 

Q2 = I + ^K( qo )(2E( qo ) - (2 - qo)K(qo)) X' (56) 

+ JLk(</ ) 3 (4(2 - qo)E(qo) - (8 - 8q + q 2 )K(qo)) A' 2 
+ 1 K^q) 5 



13 



16tt 6 (E( 90 ) - K(<zo))(E(g ) - (1 - go)K(go)) 
x (2(16 - Wq + 7q 2 )E(q ) 3 

-(2 - g ) (48 - 48q + 7q 2 )E(qo) 2 K(qo) 

+ (96 - 192g + 114gg - 18gg + ^)E(g )K(g ) 2 

-(1 - e? )(2 - gb)(16 - 16g + 9o)K(go) 3 J A 
+0(A /4 ). 

Comparing the string theory result eq. (J5fi|) to the gauge theory result eq. ljoTj) . we see 
that the 0(X' 2 ) two-loop terms precisely agree! Our joy about this amazing result would 
be even greater if this was also the case for the 0(X' 3 ) three-loop contribution. However, 
here we notice, despite a close structural resemblance (in particular, the denominators 
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of the rather involved formulas agree), that the detailed form of the gauge and string 
theory results are different. Furthermore, these analytic findings explain the numerical 
results of jTHj for the ground state gauge energies of the folded case at half filling a = \. 
There the thermodynamic ground state energies of the gauge operators were estimated 
by extrapolation of the numerically exact diagonalization of the dilatation operator for 
states up to J = 16. The method showed excellent gauge-string theory agreement (to 
about 1%) at one and two-loops, but a curious 17% deviation at three loops. Having 
at our hands the exact Inozemtsev-Bethe results eq.([51|) we can find 11 the numerically 
exact values of the three- loop ground state energies at half-filling 12 : 

Q 2 = 1 + 0.356016A' - 0.212347A' 2 + 0.212147A' 3 , (57) 

while the string result eqs. (|H4*)) . (j57)|) gives 

Q 2 = 1 + 0.356016A' - 0.212347A' 2 + 0.181673A' 3 , (58) 

We see that this explains the numerical results in Table 1 of 19J. It shows that the 17% 
deviation noticed there is not due to numerical inaccuracy of the extrapolation method 
- weak coupling gauge and strong coupling string theory results really begin to differ at 
three loops. 

The attentive reader might wonder whether the three-loop disagreement might dis- 
appear when the above constants eqs.([4i p .([44 p are chosen in a different way: Leaving 
aside the Inozemtsev chain, the general form eqs. ([40 )1 . ([43 )1 of "perturbing" the Bethe 
equations of the one-loop Heisenberg model is of course very natural. However, this is 
not the case. One checks (1) that requiring agreement between gauge and string theory 
at two- loops fixes c±, e% to precisely the values derived from the Inozemtsev-Bethe ansatz, 
that furthermore (2) no choice of C2, e2 enables one to match the solution of the deformed 
Bethe equations and the string prediction, and that finally (3) the three-loop gauge result 
eq.(jSZj) is very sensitive to the precise values of C2,e2 and the beautiful agreement with 
table 1 of [H], cf. eqs (|5Tj) . (|5l^) is only recovered if we chose the "Inozemtsev values", 
i.e. eqs. fl4"T]) . (pB|) . 

On might wonder whether the Inozemtsev chain reproduces the string result in the 
non-perturbative BMN regime where A —>■ oo in conjunction with J — > oo, while keeping 
A' = A/ J 2 small. We discussed in section 2.2 that in this case a finite scaled rapidity 
exists, cf. eq. ()14)) . Upon inversion we find 

- p(<p) = " __ log v . ( •~»9 ) 

In this context the just mentioned finding (2) (no choice of C2, ei enables one to match the 
three-loop string and spin chain results) is unfortunate since it means that the hyperbolic 
(i.e. wrappings are neglected) Inozemtsev chain at strong coupling also fails to reproduce 

11 This nice test of the validity of our method was first suggested to us by N. Beisert. 
12 It is technically convenient to use the inverse Gauss-Landen transform on the gauge result ea. (|5T)l : 
a = i corresponds to the modulus to = 0.8261148, where 2E(io) = K(£q). 
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the string prediction. Expanding eq.(59) in A' yields at three loops a potential of the same 
form as found on the right side of eq.(|4U|). One easily checks the quite interesting and 
non-trivial fact that at two loops the strongly coupled hyperbolic chain again reproduces 
the same potential as the weakly coupled chain. 

For further confirmation of these findings, let us now turn to the case of the circular 
string. 



3.4 Circular case 

Here we will be rather brief, as the details are easily filled in by combining the results and 
notations of ^Sl El E] and the methodology of the last section. The Inozemtsev-Bethe 
equation reads in this case 

'' , / / a V 2 1 , f + c A' A' 2 
dyvw)— 2 ^ = --<^log ci— + c 2 — , 60 

where we have rotated the variable (p — > tip, while the three-loop gauge energy Q 2 is 
found from 13 

Q 2 = 1 + ^5 Q 2 - ex A' 2 Q 4 + e 2 A /3 Q 6 , (61) 

with the constants Ci,c 2 ,ei,e 2 again given by eqs.fl"T ]) ..(j4"4" j) . The Q 2 k are the moments 
of the distribution density cr(</?) and obtained by 

They are generated by the resolvent 

oo „d 12 

H(<P) = ^-J2 & k ^ 2k l - e - H( V ) = 2c + <p\og^+ d<p'a(<p') / 2 . (63) 
2 j-i c + ip J c ip /2 -if 2 



k=l 

The normalization is 



ri a 
2c + I d(pa(x) = -, (64) 



where the filling fraction is a = J 2 /J > 0. The solution involves again the elliptic 
integral of the third kind, II, with modulus r = c 2 /d 2 and reads 



(65) 



1 2 / 1 , c 2 + rf 2 1 

ci A — - c 2 A 



y? 2 z V^ 4 2c 2 ^ 2 <f 2 , 
while the conditions fixing the interval boundaries c, d = c/ \fr and the modulus r are 



1 . . w 1 + r 1 w9 3 + 2r + 3r 2 1 
- i + 2K(r)c + c 1 V — - - c 2 A - = 0, (66) 

13 The sign changes in eas. (ISUll . (|^T|) w.r.t. their folded analogs eas.l)30 )l .l)35 )l are due to the rotation 
(p — ► iip, and a slightly different definition of the moments in both cases. See also |21|. 
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and 



a 



[(v9 - 1) K(r) + E(r)] + Cl A' v 2 V 7 - c 2 A' : 



/f LV ' w WJ c 2 4c 4 

(67) 

Expanding the modulus r = r + r 1 A / + r 2 A /2 and the interval boundary c, and eliminating 
all constants except ro, we find the relation 

«=±-^= + ^=§H (68) 
2 2 v ¥JK(r„) ( ' 



The three-loop gauge energy is then 

Q 2 = 1 + ^K(r ) (2E(r ) - (1 - r )K(r )) A' (69) 



2 

7T' 

16 K(r ) 5 



— K(r ) 3 (4(1 + ro)E(ro) - (1 - r )(3 + r )K(r )) A' 2 



vr6 E(r ) (E(r ) - (1 - r )K(r )) 

x ((3 + 2r + 3r 2 )E(r ) 3 - (1 - r )(8 + 3r + r 2 )E(r ) 2 K(r ) 

+ (1 - r ) 2 (7 + r )E(r )K(r ) 2 - 2(1 - r ) 3 K(r ) 3 ) A' 3 

Now we will again compare this to the string prediction, eqs. (|27j) . (|28j) . Written in 
Lagrange-inversion form, we have 



A' 



( K(t )-E(t ) V _ / E(t )-(l-t )K(t ) Y 

^ K(*)-E(t) ; V - (i - t)K(t) ; 



(70) 



tK(t) /l / K(t )-E(t ) \ 2 _ 1-t / E(to)-(l-t )K(t ) A 2 
^ 2 t K(to)V * V K(t)-E(t) J t V E(t)-(l-t)K(t) J ' 

where now the string modulus to is parametrizing the filling fraction through 

As before, power series inversion and substitution gives, from eqs.(J7UJ), 

0,2 = 1 + -^E(t )K(t ) A' (72) 



vr 2 ' 

-4K(t ) 3 ((2 - to)E(to) - (1 - *o)K(t„)) A' 2 



K(t ) e 



7r 6 E(to) 2_ (1 _ to)K(to) 2 

x ((7 - 7t + t 2 )E(t f - 7(2 - i )(l - t )E(to) 2 K(t ) 

+ (1 - i )(9 - 9t + tg)E(t )K(t ) 2 - (2 - t )(l - t ) 2 K(t ) 3 ) A' 3 

-0(A' 4 ) 
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After the Gauss-Landen transformation 



t = — v " (73) 



the string energy up to three-loop order becomes 

Q 2 = 1 + ^K(ro) (2E(r„) - (1 - r )K(r )) A' (74) 

-^K(ro) 3 (4(1 + r )E(ro) - (1 - r )(3 + r )K(r )) A' 2 
+ 4 K(r ) 5 



vr6 E(r ) (E(r ) - (1 - r )K(r )) 

x (2(7 + 2r + 7r 2 )E(r ) 3 - (1 - r )(35 + 6r + 7r 2 )E(r ) 2 K(r ) 

+ (1 - r ) 2 (29 + 2r + r 2 )E(r )K(r ) 2 - 8(1 - r ) 3 K(r ) 3 ) A' 3 

+0(A' 4 ) 

The result of the comparison between the string theory prediction eq. (|74j) and the gauge 
theory computation eq. (|69|) leads to the same conclusion as the previously treated case 
of the folded string. There is a beautiful and non-trivial agreement at the two-loop order; 
on the other hand, the three-loop terms, despite great structural similarity, are definitely 
different. Again, varying the constants c 2 ,e 2 in eqs.()60 |) .(ffiT |) does not help. 

It is instructive to consider the special case of half-filling a = 1/2. Here the solution 
becomes algebraic, and the string result reads, to all orders 

12 \/3 



A' A' 2 A 



Q 2 = v/fTA 7 = i + ___ + _ + ... (75) 

The three- loop gauge theory result, on the other hand, under the assumption of the 
validity of three-loop perturbative BMN scaling, yields 

Q 2 = l + y-^ + 0A' 3 (76) 
The three-loop contribution to the "circular" state vanishes] 

3.5 A curious observation 

It is clearly of some interest to pin down the difference between three-loop perturbation 
theory, and the string prediction. For the folded case we have 

Q2 ~ Q2 = - T |^K(g ) 5 (2E(gb) - (2 - q )K(q )) A' 3 , (77) 



and in the circular case 



Q 2 -Q 2 = - 4(1 6 ro) K(r ) 5 (2E(r„) - (1 - r )K(g )) A' 3 . (7* 

7T° 
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Curiously, this looks like a non-linear (and therefore, from the point of view of the 
perturbative spin chain, non-local) contribution to the energy. In both cases the differ- 
ence is proportional to the product of the one-loop second moment and its leading 
correction Q<p: 

Q2 ~ Q2 = -^7 Q? Qf where Q 2 = + X'Q^ + 0{\' 2 ) (79) 

But the significance or wider validity of this observation remains unclear to us at the 
moment. It would be interesting to see whether the same "recipe" also allows to account 
for the three- loop near-plane wave discrepancy of Callan et al. |28j . 
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